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The concept of an asymptotic cone of acceptance of a cosmic ray 
neutron monitor has been widely used during the past few years for the 
study of spatial anisotropies in cosmic radiation in the interplanetary 
space [ Rose and Lapointe » I96I; McCracken , 1962; Rao , et al., 

Lockwood and Razdan » 19^3 ] • This cone is defined as the solid angle 
containing the asymptotic directions of approach of cosmic ray particles 
outside the influence of the geomagnetic field that significantly contrib- 
ute to the coxmting of a ground detector [ McCracken, 1959]* .Knowing 
these cones at various stations, one can relate the Intensity obsejrved 
at a particular time to a direction in space outside the magnetosphere . 

So far, these cones have been calculated for the quiet geomagnetic field 
conditions for a single dipole model of the earth [ Lapointe and Rose , 
1961] and for a high degree simulation of the geomagnetic field 

[ McCracken. I 962 and McCracken , et ad., 1962] . For neutron monitors at 
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particular stations, differences of about 10 ° have been obtained in the 
effective directions of arrival in the letter case. However, during a 
period of geomagnetic disturbance, the so _ar wind appreciably deforms 
the geomagnetic field, confining it in a roughly ellipsoidal cavity 
with the magnetospheric boundary close to the earth on the sunward side 
and extended far out on the night side [ Sonett , et al. , 196OJ Cahil 1 
and Amazeen, 1963; Heppner , et al. , 1963]. The electromagnetic state 
of the interplanetary medium undergoes large changes at such time, 
leading to various types of observed variations in the cosmic ray inten~ 
sity, such as Forbush decreases, enhanced diurnal variation, solar 
flare Increases, etc. For the study of these variations and, conse- 
quently, the electromagnetic state of the interplanetary medium, it is 
important to know any change that might occur in the asymptotic cones 
of acceptance as a result of the distortion in the geomagnetic field. 

Thus cosmic ray trajectories were calculated nmerlcally on an IBM 709^ 
computer, and the asymptotic cones for a distorted geomagnetic field 

I 

were compared with those for a dipole field. The trajectory behavior 
near the geomagnetic cutoff for a dipole and the distorted geomagnetic 
field were also compared. 

A deformed magnetosphere was represented by Hones * model [19^3] 
where the distorted field is created by an image dipole of moment 28 
times that of earth^s equivalent dipole, placed at a distance of 28 earth 
radii along the noon meridian. The scalar potential due to these two 
parallel but unequal dipoles is given by 



28r 


U = M cos 0 + 


[(28Ee)^ + - 2(28Rg)r sin 0 cos 


-J 


( 1 ) 


where 


M moment of earth’s dipole (8.1X10^® gauss cm®) 

0 colatitude 

9 azimuthal angle measured aroiuid earth’s dipole from the noon 

meridian 

r radial distance from the center of the earth 

Eg mean radius of the earth 

This leads to a magnetic field configuration in which the earth’s field 
is completely confined in an approximately ellipsoidal cavity extending 
7*56 Re "the solar direction and 23.91 Re in the night direction. 

A cross-sectional view of the deformed magnetosphere in the meridian 
plane is shown in Figure 1. It is found that an ellipse with major 
and minor aoces of 15*736 Rg nnd 14.2^0 Rg fits the meridian plane 
trace fairly well. By then making a simplifying assumption that the 
boxmdary of the deformed magnetosphere is an ellipsoid of revolution, 
we arrive at the boundary equation for the ellipsoidal magnetosphere 
given by 



(r sin 0 cos cp + 8.176 Re)^ (r cos 0) (r sin 0 sin cp)^ 

+ _ — ^ ^ ^ ^ ^ ^ _ 


(2) 


(15^736 Re)^ (li )-.24 Rg)2 (ik.2l|- Rg)® 

It is clear from equation (l) that, by changing the strength of the image 
dipole and its distance from the center of the earth, we can obtain 
other types of ellipsoidal models for the magnetosphere where the 
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position of the "boundary on the solar side and the ni^t side would he 
different . However , the above model which would approximately fit the 
experimental obseirvations during a severe geomagnetic distvirbance would 
be reasonable for the present calculations . , 

The asymptotic directions of cosmic ray particles outside the 

'A 

magnetosphere and arriving at a particular geomagnetic latitude were 
calculated by considering an equivalent problem of the motion of negative 
particles of the same rigidities moving in the reverse direction. The 
equation of motion 



where the symbols have their usual meaning, was solved in a spherical 
coordinate system by the Runge-Kutta-Blum method [ Blum, 1957] • When the 
particles reached the boundary of the distorted magnetosphere, as given 
by equation (2), the trajectory direction was noted in terms of asymptotic 
geomagnetic latitude (A) and asymptotic geomagnetic longitude (t)^~^as- 
ured eastward from fetation meridian. Asymptotic directions were calcu- 
lated for various particle rigidities for the dipole field, and the 
distorted field at station azimuth angles of 0°, 90° l80°, and 270°. 

For the dipole field the integration of the equation of motion was ter- 
minated at 25 earth radii, beyond -Hhich the geomagnetic deflection siiffered 

i 

by a particle is negligible . | 

In Table la, A,i|f values are presented for Vertically incident 
particles near the cutoff at 50° geomagnetic latitude. Particle 
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riglditlee from 2.7 to 3*2 bv at steps of 0.1 bv are considered which 

mostly cover the penxjmbral region at this geomagnetic latitude (Stormer | 

cutoff at 50° = 2.5^ bv). In Table lb, the A,\lr values are given ^Tables la | 

Xauid lb S 

only for a 3 bv particle arriving at 50° geomagnetic latitude along ^S, f 

the vertical and zenith angles of l6° and 32° in the geomagnetic North, ! 

South, East, and West planes. A point which may be noted in these tables ■ 

t 

is that some particle trajectories in a distorted field, which are allowed | 

t 

• in one longitiide ( i .e • , <p = l80° ) , may be prevented from entering in | 

I 

another. QThis indicates that during a geomagnetic disturbance there jl 

might occur changes in the alloved and forbidden cones in the penumbral I 

region as a function of local time of the station. Recently, Shea , et al . | 

[1964] have studied the trajectories in penumbral region at various geo- | 

magnetic latitxides at intervals of 0.01 bv, using a high degree simulation | 

of the quiet geomagnetic field. For every forbidden trajectory, they 

have increased the StSrmer cutoff at a particular station by 0.01 bv 

and arrived at what .they call the effective geomagnetic cutoffs which give i 

a better fit to the e3q>erimentally observed variations in the cosmic ray 

intensity as- a function of geomagnetic latitude and longitude. From the 

results presented in ITables la and lb, it appears then that during a 

geomagnetic disturbahee these effective cutoffs would vary along different 

longitudinal planes, which would cause local time effects in especially 

low energy cosmic ray intensity observed at balloon heists or in low 

i 

altitude satellites. No attempt has been made here to calculate these 
effective cutoffs, since for that purpose it would be appropriate to 

I 



consider an aocxirate model of the distorted geoma^etic field where the 
higher harmonics of the surface field are also taken into account. It is 
also apparent from. Tables la and lb that except for the midnight meridian ^ 
there are more nonentrant particles in the distorted than in the dipole 
geomagnetic field. This result, though derived from a limited number of 
trajectories, is in agreement with the conclusion drawn by Akasof u , 
et al . [1963] that the limitation of the radial extent of the geomagnetic 
field due to the solar wind cannot produce reduction in the geomagnetic 
cutoffs . 

To determine asymptotic acceptance cones, the asymptotic directions 
for particles of various rigidities were calculated for arrival at eeirth 
from the vertical and from the geomagnetic North, South, East, and West at 
zenith angles of 16® and 32°. The rigidity intervals from geomagnetic 
cutoff to an upper limit of 200 bv were chosen such that there were 
relatively small changes in asymptotic directions from one rigidity to 
the next. This is not strictly true near the station cutoff as discussed 
later. The asyn^it'otic directions corresponding to each rigidity band and 
the arrival cone were -given proper weighting for the cosmic ray primary 
rigidity spectrum taking into account the yield functions of neutrons at 
various atmospheric depths [ Webber and Quenby , 1959 3 ^ and the known zenithal 
response of neutron monitors. For further details of the method of calcu- 
lations, see Lapointe and Rose [1^9^13, McCracken [19^23, and Rao, et al . 
[19633* The asymptotic cones of acceptance of* sea level neutron monitors 
have been calculated in a dipole and the distorted geomagnetic fields at 
three representative geomagnetic latitudes, 0°, 50° ^ and 70°, for azimuthal 
angles of 0°, 90°, l80°, and 270°. The results are presented in 
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• Table 2 where an asymptotic cone of acceptance of a neutron monitor is ^able 2 
represented by A and 'if, the weighted mean values of asymptotic geo- 
magnetic latitudes and longitudes of particle trajectories, and 

and their root -mean squares, which give an idea of the width of 
the acceptance cone. A niaxiraum difference of about 6® is observed in 
these values at various geomagnetic latitudes when the acceptance cones 
in a dipole field axe compared with those calculated in a distorted 
geo m ag n etic field at various azimuthal angles. These differences do not 
seem to vary in a definite manner with the geomagnetic latitude or longi- 
tude. We believe these axe largely due to errors introduced in the 
calculation, particularly because of trajectories neax the station cut- 
off . The asymptotic directions of particles change very rapidly there 
with as small ,a change in rigidity and the arrival direction as 0.1 bv / 
and 1°, resi)eotively . It is therefore difficult to choose a paxtlcle / 
trajectory which would represent the average behavior of the whole 
rigidity band near the cutoff for various cones of arrival, tinless 
rigidity bands, and ‘arrival cone widths of less than 0.1 bv and 1®, 
respectively, are considered. Because of the laxge amount of time consumed 
on the computer such fine intervals have not been considered here. Besides, 
the errors involved ^ the specific yield functions of cosmic ray particles 
and in our knowledge of the zenithal response of neutron monitors would 
also cause uncertaiirties in the calculated asymptotic cones of acceptance 
which would easily be of the order of It may therefore be concluded 

that during magnetically distvirbed periods, the change in the geomagnetic 
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field configuration due to solar wind does not by itself bring about any 
appreciable change in the asymptotic cones of acceptance of cosmic ray 
neutron monitors at sea-level stations beyond the possible errors involved 
in the method of calculations. A model calculation for a high altitude 
station (atmospheric depth 680 g/cm^), at the geomagnetic eqviator, shows 
that thOTJgh there is a shift in the values, the change from dipole 
to distortion field is negligible. However, dui’ing magnetically disturbed 
periods, there may be other effects such as a ring cxorrent around the 
earth, a change in the cosmic ray rigidity spectrum, or a cheinge in sta- 
tion cutoffs, which might modify these acceptance cones. The former . 
effects have been discussed by McCracken [I962], McCracken, et [1962], 
and Ran, et ad. [I963l* The latter effect would again be negligible for 
small changes in station cutoffs for at least those neutron monitors 
below the atmospheric depth of 680 g/cm^. 

In the present calculations a highlj’’ distorted geomagnetic field was 
considered where one would expect to find maximum changes in the a^mp- 
totic cones of acceptance at various latitudes and longitudes. The nega- 
tive results therefore imply that there would be no significant effects 
for less severely distorted geomagnetic fields. Further, a steep cosmic 
ray rigidity spectrum corresponding to a minimum solar active period of 
195^ [ Webber and Quenby , 19591 was used for giving proper weighting to 

i 

asymptotic directions of various particle rigidities. This implies that 
for cosmic ray variations with flatter spectrums, such as Forbush decreases 
and the diurnal variation, the change in the asymptotic cones of acceptance 
due to geomagnetic field distortion would again be insignificant. However, 


for solar flaxe type variations in cosmic rays, asyir^rtotic acceptance 
cones may undergo changes due to geomagnetic field distortion, since 
the spectrum for such variations is much steeper, thereby leading to 
a greater bias for asymptotic directions of low rigidity particles 
which undergo large changes. A model spectrum of the type R~® 

(R = Rigidity) was used to calculate asymptotic cones of acceptance 
of neutron monitors at 680 g/cm^ atmospheric depth at geomagnetic lati- 
tudes of 50° and 70° for the dipole and the distorted geomagnetic field 
along the noon meridian. The weighted mean values of asynptotic geo- 
magnetic latitudes and longitudes presented in Table 3 show that while 
there are large changes in t at 7 ^ = 50° the changes at = 70 ° 
occur mainly in A. It should be recognized though that the errors 
associated with if and A are also larger here because of more bias 
toward asymptotic .directions of particles close to the geomagnetic cut- 
off . Extending these results to detectors at sea level and those at 
the top of the atmosphere, one would expect smaller changes in asymptotic 
cones of acceptance in the former case but much larger changes in the 

I 

latter case. It therefore seems that when anisotropic effects of cosmic 
ray solar flare increases are studied by means of balloon borne detectors 
or low altitude satellites, it is important to take into account the 
geomagnetic field distortion in calculating asymptotic cones of accept- 
ance . For each Individual event , an appropriate rigidity spectrum and 
the geomagnetic field distortion would have to be considered. Besides, 
it would also be important to take into accoimt the higher harmonics of 
the surface geomagnetic field which, together with solar wind distortion 
effects, may cause a larger change in the asymptotic acceptance cones at 



certain stations 
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TABLE 1 .- Asymptotic geomagnetic latitudes and longitudes of 
particles in dipole and disto?rxed field 
(a) Vertically arriving cosmic rays near cutoff at 50° 
geomagnetic latitude 


Rigidity 

Dipole field 


Distorted 

field 


Cp = 0 ° 

90° 

180° 

270° 

2.7 "bv 

A 92.3° 

i 304 . 0 ° 

No entry 

95.6° 

4379.1° 

92.8° 

525.8° 

No entry 

2.8 bv 

^ No entry 

93.1° 

918.40 

No entry 

94.7° 

821.7° 

No entry 

2.9 bv 

A 89.0° 

t 353.90 

100.2° 
656 .5° 

80.0° 
489 .6° 

84 . 6 ° 
383 .2° 

99.2° 
895 .9 

3-0 bv 

A 101 . 9O 
Ilf 254.6° 

No entry 

No entry 

90.8° 
383 .5° 

No entry 

3.1 bv 

A 73 . 1 ° 

t 179.1° 

102.3° 

445.1° 

102.8° 

426 . 0 ° 

91.1° 

379.1° 

86.8° 

599.0° 

3*2 bv 

A 81.7° 
\jr 151.8° 

82.3° 

229.7° 

74.1° 

177.0° 

78.6° 

152.1° 

81.0°^ 

509.7- 


t 
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table 1 .- Asymptotic geomagnetic latitudes and longitudes of 
particles in dipole and distorted field - Concluded 
(b) A 3 Tyv particle arriving at 50 ° geomagnetic latitude along 
the vertical and zenith angles of l 6 ° and 32 ° in the 
geomagnetic North, South, East, and West planes 


T,. Arrival 

Blgidlty direction 

Dipole 


Distorted field 


field 

Cp = 0 ° 

90° 

180° 

270° 

Vertical 

A 101 . 9 ° 

t 254.6° 

No entry 

No entry 

90.8° 

383.5° 

No entry 

l 6 ° N 
32 ° N 

A 99.4° 
\lf 402 . 6 ° 
A 76.3° 
^ 168.6° 

No entry 

91.9° 

267.5° 

No entry 

93.2° 
234 . 4 ° 

83.5° 

271.6° 

97.9° 

275.5° 

No entry 

103.5° 

382 . 1 ° 

3 hv l 6 ° E 

>• 32 ° E 

A 88 . 4 ° 

t 1092.7° 

^ No entry 

90.4° 
274.5° , 

No entry 

89.0° 
462 . 2 ° 
^ 87.7° 
620.2° 

No entry 

95.5° 
445 .8 

98 . 8 ° 

387.0° 

No entry 

l 6 ° S 
32 ° S 

A 78.7° 
t 209 .0° 

A 75 . 2 ° 

t 193 -2° 

No entry 
No entry 

87 . 6 ° 

408 . 8 ° 

100 . 1 ° 

528 . 5 '. 

88.5° 
362 .2° 
92.0° 
290.8° 

No entry 

96.3° 

635.2° 

l 6 ° W 
32 °. W ' 

A 73.3° 

V 174.0° 
A 77.6° 

V 156.0° 

79.6° 

721.0° 

108.1° 

287.9° 

87.1° 

564 . 0 ° 

72.5° 

190.8° 

87.6° - 
244 . 0 ° 
71.6° 
167.2° 

76.7? 

346 . 2 ° 
106.8° 
293 .1° 
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TABLE 2 .- Mean asymptotic .geomagnetic latitudes and longitixdes 
their root -mean squares ( ^ asymptotic 
cones of acceptance of cosmic ray neutron monitors at 
various geomagnetic latitudes and longitudes^ for the 
cases of dipole and the distorted geomagnetic fields 


Geomagnetic 

latitude 

Dipole 

field 


Distorted 

field 


cp = 0° 

90° 

180° 

270° 

0° 

A 

0° 

0.1° 

0° 

0° 

0° 



70.2° 

66.1° 

70.2° 

67.4° 

63.8° 



10.0° 

9.8° 

10.5° 

10 .3° 

10.5° 



50 . 5 ° 

53 - 1 ° 

50.2° 

47.8° 

51 .8° 

50 ° 

A 

6 . 9 ° 

11 . 8 ° 

7.4° 

7.1° 

9.9° 


t 

49.3° 

43.3° 

49.6° 

49.6° 

42 . 5 ° 


°A 

22.3° 

21.5° 

22.3° 

21 , 8 ° 

26.1° 



31.2° 

34.6° 

31.4° 

36.5° 

36.7° 

70 ° 

A 

{ 

47.0° 

50.4° 

48 . 6 ° 

46 . 9 ° 

48 . 6 ° 



24 . 3 ° 

21 . 2 ° 

25 - 1 ° 

26.0° 

22 . 4 ° 



11 . 6 ° 

10.5° 

11.3° 

12 . 2 ° 

11.3° 


% 

18 . 4 ° 

19.0° 

18.7° 

18.4° 

18.3° 



r 



Figure Title 


Fig. 1.- Meridian across section of deformed magnetosphere [Hones ^ 1963]* 




